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ABSTRACT 



Dust-grain growth and settling are the first steps towards planet formation. An 
understanding of dust physics is therefore integral to a complete theory of the 
planet formation process. In this paper, we explore the possibility of using the 
dust evaporation front in YSO disks ('the inner rim') as a probe of the dust 
physics operating in circumstellar disks. The geometry of the rim depends sensi- 
tively on the composition and spatial distribution of dust. Using radiative trans- 
fer and hydrostatic equilibrium calculations we demonstrate that dust growth 
and settling can curve the evaporation front dramatically (from a cylindrical ra- 
dius of about 0.5 AU in the disk mid-plane to 1.2 AU in the disk upper layers 
for an AO star). We compute synthetic images and interferometric visibilities for 
our representative rim models and show that the current generation of near-IR 
long-baseline interferometers (VLTI, CHARA) can strongly constrain the dust 
properties of circumstellar disks, shedding light on the relatively poorly under- 
stood processes of grain growth, settling and turbulent mixing. 

Subject headings: young stellar objects — circumstellar disks — radiative transfer 
— Monte Carlo codes — dust sublimation — grain evolution — interferometry 



1. Introduction 



Advances in long-baseline near-infrared interferometry have made it possible to study 



the inner regions of circumstell a r disks at sub-AU scales. Early results (iMillan-Gabet et al. 



19991 . I2OOII : iTuthill et all I2OOII : iMonnier fc Millan-Gabetl 120021 1 have shown that the dust 



disk gets truncated at a finite radius (determined by the luminosity of the central star 
and dust sublimation temperature) with in which the tenaperature is too high fo r dust to 
survive. The truncated disk forms a 'rim' ( Natta et al.ll200ll : iDuUemond et al.ll200ll . hereafter 
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DDNOl) that subtends a significant solid-angle around the star and intercepts stellar photons, 
re-radiating predominantly in the near-IR. This rim naturally explains the near-IR excess 
observed in Herbig Ae systems. 

Despite its success in accounting for the near-IR excess, the vertical rim model intro- 
duced in DDNOl has had a few weaknesses from the very outset. The flux received from a 
vertical rim is highly viewing-angle dependent, since the projected line-of-sight area of the 
rim changes rapidly with inclination angle. This implies that the near-IR spectral energy 
distributions (SEDs) of Herbig Ae/Be stars should be a strong function of inclination angle, 
but observations show that most Herbig Ae stars ha.ye similar near-IR ex cesses, independent 
of their inferred inclination angles (iNatta et al.ll200ll : lDominik et al.ll2003l ). Furthermore, the 
surface brightness distribution of a vertical rim will become highly asymmetric at inclination 
angles different from face-on. This asymmetry sh ould be detectable by interferometers in 
the form of a closure-ph ase signal (IMonnieii 12000 . and references therein), but recent stud- 
ies (IMonnier et al.ll2006l ) do not detect significant closure-phase signals indicating that the 
geometry of the inner rim is not a vertical wall. 



It appears that a mechanism that curves the inner rim is needed. Ilsella fc Nattal (2005) 
(hereafter IN05) studied the impact of density-dependent dust sublimation temperatures 
on the structure of the rim. Close to the m id-plane, where the densities are high, the 
dust sublimation temperatures are also high (jPoUack et al.l Il994l ). while as the densities 
begin to fall rapidly in the vertical direction the sublimation temperatures drop. The IN05 
calculations show that this causes the inner rim to curve away from the surface of the star. 
Isella et al.l (120061 ) found that a curved-rim model provides a substantially better fit to the 
visibilities and near-IR SEDs of a few Herbig Ae systems. 

In this paper we look at another mechanism that could curve the rim on larger scales than 
the ones produced by a density- dependent dust sublimation model — namely, dust growth 
and settling. In the absence of strong turbulence, dust grains will settle towards the disk 
mid-plane with large grains settling lower and having a smaller scale height than small 
grains. Dust grains will also grow rapidly in the mid-plane due to the m uch higher densities 
(IDubruUe et al.l Il995l : iDuUemond fc DominikI |2004J : iTanaka et al.l |2005| ) . The opacity in the 
mid-plane is thus dominated by large grains whereas small grains dominate the opacity higher 
in the disk photosphere. This transition in the nature of the opacity source along with dust 
sublimation curves the inner rim as described in Section 4. Near-IR interferometers with 
baselines several hundred meters long (e.g. CHARA, VLTI) can distinguish between the 
model presented in this paper and the IN05 model, shedding light on dust physics operating 
at sub-AU scales. 



We have computed the structure of the dust evaporation front in a self-consistent manner 
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for the first time in a 2-D geometry (photon scattering has not been in cluded; for effects of 
scattering on the disk structure and SEP see iDuUemond &: Nattall2003l ) using the radiative 
transfer code TORUS (iHarried |2000| ) . In the following section we describe the radiative 
transfer code, while in section 3 we perform a consistency check by comparing the results 
for a TORUS generated density-dependent dust sublimation model with IN05 results. In 
Section 4 we present results for the dust segregation model, followed by a summary of the 
results and conclusions in Section 5. 



2. The Monte Carlo Radiative Transfer Code - TORUS 



2.1. Description of grid and disk structure calculations 

The calculations in this paper were performed using the TORUS Monte-Carlo radiative- 



transfer code (lHarriesll200C 



3 paper were periormed using tne torus Monte-Uario raaiative- 
Harries et all2004 : Kurosawa et al. I l2004h . Radiative equilibrium 



is computed using Lucy's (ILucy Ill999l ) algorithm on a two-dimensional, cylindrical-cartesian 
adaptive-mesh grid. Storing the opacity and temperature information on an adaptive mesh 
has particular advantages for the problem considered here, since the accurate determination 
of temperatures and subsequent SEDs requires an adequate resolution of the effective disk 
photosphere, whose position is changing spatially both as the dust is sublimated and as the 
vertical structure of the disk changes as it is iterated towards hydrostatic equilibrium. Using 
adaptive-mesh refinement (AMR) we are able to subdivide cells near the disk photosphere 
as the calculation proceeds, ensuring the resolution is maintained at each step. 

The initial density structur e for the disk calcula t ions i s based on the canonical descrip- 
tion of the a-disk developed by IShakura fc SunyaevI (119731 ) . viz 



p(r, z) 
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where r is the radial distance in the mid-plane, tq is some characteristic radius, z is the 
distance perpendicular to the mid-plane, and h{r) is the scaleheight, given by 



h{r) 



hn I - 



(2) 



with parameters of a = 2.625 and f3 = 1.125, giving a radial density dependence of the 
surface density of S(r) cx r^^-^. The AMR mesh is divided such that the cells are loga- 
rithmically spaced in the radial direction and that there are at least 7 cells in the vertical 
direction per disk scale height h{r). A sweep of the grid is then made to split cells further 
around the tssoo = 1 (T5500 is optical depth calculated at 5500 A) surface using the following 
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algorithm: for each pair of adjacent AMR cells the individual optical depths across the cells 
are calculated (ri and T2 say, with ti referring to the more optically thick of the two). The 
more optically thick cell is subdivided providing that T2 < Tmin and ti > r^ax, where r^i^i 
and Tmax are parameters (typically 0.1 and 0.5 respectively). Should any cell in the grid be 
subdivided, the sweep is repeated until no new cells are added to the mesh. 

Once the temperature (we assume that the disk is in local thermodynamic equilibrium 
and gas and dust are thermally coupled) and dust sublimation (see next section) struc- 
tures have converged using the Lucy algorithm, the vertical disk structure is modified via 
the equation of hyd rostatic equilibrium following a similar algorithm to that detailed by 



Walker et al.l (120041 ). Briefly, and assuming disk mass is negligible compared to the central 



star, the equation of hydrostatic equilibrium is 

^ = -P3^ (3) 

where P is the pressure and Qz is the vertical component of the stellar gravity. Adopting an 
ideal gas equation of state P = pkT/ fimn (with /i = 2.3 for a standard molecular hydrogen- 
helium mixture) we get 

dp 1 If gzumn ^ dr\ 



dz p T \ k dz J 

The above equation can be solved for p{z) numerically on the AMR mesh since the verti- 
cal temperature structure is known from the radiative-equilibrium calculation. The vertical 
density structure is then renormalized to conserve the radial surface density profile detailed 
above. A self-consistent calculation for dust sublimation (see §2.2p and disk temperature fol- 
lowed by a hydrostatic equilibrium calculation is repeated until the disk density structure has 
converged. Convergence is typically achieved in four iterations. Images and SEDs are sub- 
sequently calc ulated using a separate Monte Carlo algorithm based on the dust emissivities 
and opacities ( Harrieslboool ). 



2.2. Implementation of Dust Sublimation. 

In the models discussed here the shape of the inner rim is controlled by dust sublimation. 
Typical mid-plane densities in the dust sublimation region are of the order of 10~^ — 10^^° 
gcm~^ in Herbig disks. In the mid-plane, an optical depth of unity for visible light photons 
(measured radially inward into the disk from the rim edge) corresponds to a length scale of 
~ 10^^ — 10^^ AU. Therefore if we started sublimation runs with a grid that had the full 
dust opacity present, convergence towards equilibrium rim shape would be extremely slow. 
This is very similar to the problem encountered in photoionization codes, where neutral 
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hydrogen opacity to ionizing photons is so large that starting from a neutral grid would 
take an impractica l amount of CPU time for the Stromgren sphere to propagate outwards 



f lWood et al.ll2004f ) 



The solution that we adopt is also analogous to the one used in photoionization codes 
JWood et alibooi ). The dust content is first reduced to a very low value in the computational 



grid for the circumstellar disk, to make each of the grid cells optically thin. Stellar photons 
then propagate through the disk and the temperature of grid cells is determined. Dust is 
added to cells that are cooler than the sublimation temperature in small steps of r. The step 
size is a r of 10~^ (computed at 5500A) for the first five dust growth steps. The step size is 
then increased logarithmically, first to 10~^, then to 10~^ and so on until the appropriate gas 
to dust ratio is reached in each grid cell. The grid cell temperatures are recomputed after 
every dust growth step and the process is repeated until the shape of the dust sublimation 
region converges. 

Figure [1] illustrates the dust-growth process in our computational scheme. The color 
scheme shows integrated optical depth, measured along lines perpendicular to the disk mid- 
plane. Panel-(a) is a grid for which the optical depth across each cell is 10^^. Panel-(b) 
is a cross section of the grid after the first dust growth step where r across each cell is 
now 10^'^. The cells have a r of 10~^ across them in panel-(c). The optical depth across 
each of the cells that lie away from the rim surface is 10 in panel-(d), whereas cells at the 
outer surface of the rim have been smoothed to resolve the r = 1 surface. The geometry 
of the rim can considered to be 'sublimation converged' in panel-(d). Any further addition 
of dust does not significantly change the distance of the r = 1 rim surface from the star. 
An appropriate amount of dust is put into each of the cells to reach a predetermined gas to 
dust ratio in the final dust growth steps. A hydrostatic equilibrium calculation as described 
in §2.11 is then performed. The dust growth and hydrostatic equilibrium calculations are 
repeated until convergence for the entire disk structure is reached (typically 4 iterations). In 
the next section we describe checks that ensure that this process gives accurate results. 



3. Code testing against IN05 results 



IN05 calculated inner rim models with astronomical silicate dust ( IWeingartner fc Draine 
200ll ). where the shape o f the rim is set by th e fact that the evaporation temperature of dust 



depends on gas density (jPoUack et al. 1994 ). In this model, silicate grains sublimate at a 
higher temperature compared to other grains and hence fix the rim location. For silicate 
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dust, the evaporation temperature Tgvp can be parameterized as 

rp _ ^ Pgas(^, z) 

L Ig cm~'^ J 

where G = 2000 K, 7 = 1.95 x 10"^ and Pgas is the density of gas in gcm~^ (see IN05 eq. 
[16]). IN05 showed that this dependence of Tevp on gas density curves the inner rim. 

In order to test the numerical scheme described in §2.21 we compared TORUS results for 
the rim with the IN05 results (kindly provided by A. Isella). Table [1] describes the properties 
of the star and the circumstellar disk used in the comparisons. Figure [2] shows the shape of 
the rim for large (1.2/im) and small (0.1/im) grains. The rim here is defined as the r = 1 
surface (for A = 5500A), calculated along radial lines from the central star towards the 
disk. Figure [3] compares the fraction of stellar luminosity re-emitted by the rim in infrared 
wavelengths (integrated over 1.25-7/im) as computed by IN05 and TORUS. 

Figures [2] and [3] show that that the TORUS results are consistent with the IN05 calcula- 
tions at better than 15% level, giving us confidence in our Monte Carlo (MC) computations. 
We note that benchmark comparisons for MC codes, studying radiative transfer in optically 
thick disks (Pascucci et al. 2004), show a similar level of agreement among the various MC 
codes. 



4. Rim curvature due to grain growth/settling 

4.1. Model Description 

The rim in the IN05 model is curved because of a density-dependent dust-evaporation 
temperature. In this section we describe an alternative model in which the dust sublimation 
temperature is kept fixed at 1400 K and rim curvature arises entirely because of dust growth 
and settling. 

Small grains are much less efficient at cooling than large grains. If only single-sized grains 
are present, then at a given distance from the star, small grains will have a higher temperature 
than large grains. Under the assumption of a good temperature coupling between gas and 
dust (at densities characteristic of the rim re gion, this is a good approximation even in the 
disk photosphere; iKamp fc DuUemondl I2004J ) . a mixture of small and large grains will be 
hotter than the case where only large grains are present. Combined with dust sublimation, 
this means that the location and shape of the inner rim will change depending on the relative 
fraction of small and large grains. 



A number of observational (ID'Alessio et al.ll200ll . I2OO6I : IChiang et al.ll200ll : iRettig et al. 
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20061 ) and theoretical (iDullemond &: Dominikll2004j : iTanaka et al.ll2005l ) results have shown 
that dust properties in circumstellar disks evolve with time. Small grains tend to coagulate 
into larger grains and settle towards the disk mid-plane, with the evolution of dust grain sizes 
occurring most rapidly in the disk mid-plane where the densities are high. In the absence 
of strong turbulence (the nature and strength of turbulence in circumstellar disks is not 
well understood), large grains will have a smaller scale-height than small grains. Thus grain 
growth and settling provide a natural means for changing the fraction of large and small 
grains as a function of height ('dust segregation') from the disk mid-plane. This will curve 
the rim as described in the previous paragraph. 

As an illustration, we look at a model (henceforth 'dust segregation model') with silicate 
dust (jWeingartner &: Drainell200ll ). In order to simulate dust settling we fix the large (1.2/im) 
grain scale height to be 60% of the gas scale height. We do not settle the small grains (O.l/im) 
and their scale height is fixed to be the same as the gas scale height. As grains grow, the 
net mass of small grains decreases, while the net mass of large grains increases: we fix the 
mass of large grains to be ~ 9 times the mass of small grains. All the the dust mass resides 
in the 1.2 & 0.1 fim grain components and total mass of dust in the disk is chosen to be 
1/lOOth of the total gas mass. The rim gets curved over a wide range of parameter choices. 
The curvature can be made stronger by increasing the relative fraction of small grains. 



The above description leads to the following expressions for dust density: 

1 



pi.2 = 1.5 X 10 Po^asexp 



2 (0.6/i)2 



(6) 



po.i = 10 Vo-gas exp 

where pi.2 is the density of 1.2yum grains, po.i is the density of 0.1pm grains, Po_gas is the gas 
density in the mid-plane, z is distance perpendicular to the disk mid-plane and h is the gas 
scale height. The star and disk parameters are chosen as in Tabled] (the parameters are the 
same as in the IN05 models of section 3). 



Iz^ 
2/^ 



(7) 



4.2. Geometry of the Dust Segregation Rim 

4.2.1. Rim Shape 

Figure [2] shows the shape of the inner rim for the dust segregation model. The curvature 
of the rim depends on the relative density fractions of large and small grains and the width 
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depends on the sizes of the two dust components. A major finding here is that the dust 
segregation model can curve the rim on much larger scales than a density-dependent dust 
sublimation temperature (IN05) model. The r = 1 surface is nearly vertical at the inner 
edge of the dust segregation model because of the fixed, density-independent dust sublimation 
temperature. 



4-2.2. Analytic Estimate for Rim Shape 

The most important dust parameter that determines the rim shape is e{z), defined as 
the ratio of Planck mean opacities at the dust sublimation temperature and the stellar effec- 
tive temperature (IN05). Using equations (6) and (7) and values of Planck mean opacities 
for 0.1/im and 1.2/im grains, e (see appendix for a more general expression) for the dust 
segregation model is given by 



0.42 + 0.015 exp 



0.89^-' 



0.76 + 0.181 exp (^) 



(8) 



Figure shows e as a function of mid-plane height. Corresponding to this e, an effective 
dust size which has the same e as the dust mixture, is plotted in Figure |Dd. 



i?evp [AU] = 0.034 



[AU] 



1500 K 



evp 



Calvet et al. 


1991 
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sella et al. 


2006) 




- 1/2 



(9) 



Te^p is the dust evaporation temperature, is luminosity of the star in solar units and B is 
a dimensionless diffuse heating term. B characterizes the optical thickness of the rim to its 
own thermal emission. 



The gas scale height h at Revp is estimated to be (see eqs. [7] and [8] from lChiang &: Goldreich 



19971 ) 



h [AU] = i?evp [AU] 



gas 



4 X lO^K 



Me 



evp 



4.64 X 10-3 AU 



nl/2 



(10) 



where Tgas is the temperature of the bulk of the gas (assumed to be isothermal) and M^, is 
mass of the star in solar masses. 

Combining equations (8), (9), and (10), substituting for stellar parameters from Tabled] 
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and setting Tevp = 1400 K, the height of the evaporation front (see appendix for a general 
result) as a function of distance from the star along the disk mid-plane becomes 



^[AU] = 1.06i?evp [AU] 



In 



5.87AU" 



'Rl^ - 0.425 



0.76 



0.181 + 0.015S -0.21AU' 



''evp 



1/2 



gas 



1 X lO^K 



R, 



cvp 



4.64 X 10 
(11) 



-3 



AU 



Figure [2] shows an ana lytic estimate for t he rim shape (we assume a gas temperature of 



1000 K and have set B = 2. ICalvet et al.lll99ll ). The analytic rim agrees with the TORUS dust 



segregation rim at the 30% level. The match at the inner edge of the rim can be improved 
by tuning the diffuse heating parameter B, but the match at the outer edge remains poor 
because the analytic estimate is not self consistent and does not take into account the 
transition between the rim and the 'shadow' (see section §4.3.2p region beyond it. 



4.3. Observables 

4.3.1. SED 

Figure [5] shows the infrared emission from the inner rim and the outer flared-disk com- 
ponent. There is significant emission over the stellar blackbody (denoted by the continuous 
line) from 1.5/im onwards. The strength of the rim emission is approximately determined 
by the temperatures along the rim at optical depths of 2/3 (calculated at the wavelength 
of emission). Figure 3 in the IN05 paper shows that at a given optical depth, the small 
grain rim has a lower temperature than the large grain rim. Because of this temperature 
difference, near-IR emission in the small grain model gains strength at longer wavelengths 
compared to the large-grain and dust-segregation models. 



4-3.2. Images and Visibilities 

In order to illustrate the differences between the IN05 and dust segregation models we 
have computed synthetic images and visibilities for the different models. At 2.2/im, all the 
emission comes from the rim and the central star. The panels in Figure show 2.2/im 
rim surface brightness for two inclination angles. As discussed previously, and is evident 
in Figure O the dust segregation rim is thicker than the IN05 rims. The dust segregation 
models also distinguish themselves observationally from the IN05 rims by their visibility 
curves. Figure [7] shows plots of 2.2/im visibilities calculated along the disk major and minor 
axes. The visibilities at long baselines are more oscillatory in the IN05 models. The IN05 
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rims are also more skewed at large incli nation angles which would give them a stronger 
closure phase signal (IMonnier et al.ll2006l ). 



At 10.7/im, emission comes from two distinct regions of the circumstellar disk as seen 
in Figure [HI The rim has a high surface brightness at 10.7/im, but only contributes about 
30-50% of the integrated light (more precisely 38% in the large grain model, 55% in the dust 
segregation model and 32% in the small grain model). In contras t to 2.2u,m. emission, a large 
fraction of the emission a t 10.7/im comes from the flared disk (IKenyon fc HartmannI 119871 : 
Chiang fc GoldreichI 119971 ) regi on on larger scales. B etween the inner rim and the flared 



disk, there is a 'shadow region' (|Dullemond et a. 



directly. As discussed in Ivan Boekel et al.l (120051 ) . the two distinct regions of emission cause 



200 ll ) which is not illuminated by the star 



a knee in the visibility curves (Figure [HI). The size and shape of the knee is affected both by 
the inner rim as well as the flared disk structure (Ivan Boekel et al.ll2005l ). 



van Boekel et al.l (120041 ) studied the silicate feature in Herbig disks using mid-infrared 



interferometry (with VLTI-MIDI), and showed that dust on AU scales in Herbig disks has 
undergone significantly more thermal processing as compared to the dust at tens of AU scale. 
The discussion in this section demonstrates that the present generation of long-baseline near- 
and mid-infrared interferometers can be used to distinguish between inner-rim models, thus 
providing information on dust processes that take place on even smaller scales (sub-AU). 



Discussion 



The DDNOl model was introduced to explain the near-IR excess observed in Herbig 
Ae systems. The near-IR excess in this model comes from a vertical inner dust rim which 
is frontally illuminated by the central star. As discussed in section [H the emission from 
a vertical rim geometry becomes highly skewed even at mod erate line of sight inc lination 
angles of the disk. Recent closure phase results on YSOs by iMonnier et al.l (120061 ) do not 
find any evidence for the l arge skewness in near-IR emission predicted by the DDNOl model. 
The IMonnier et al.l (120061 ) results favor the existence of curved inner dust rims and discount 
the vertical inner wall of DDNOl. The first and only detailed calculation existing in current 
literature for rim curvature was performed by IN05, where the rim gets c urved because of a 
density dependence of dust sublimation temperatures. Initial fits to data (llsella et al.ll2006l ) 
seem to suggest that the curved rim-model provides a better fit to the visibilities and near-IR 
SEDs of a few Herbig Ae systems, although detailed comparisons have not been made yet. 



In this paper we introduced a new mechanism based on dust growth and settling to 
curve the inner rim. As a consequence of growth and settling, opacity in the rim mid-plane 
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will be dominated by large grains whereas small grains dominate the opacity higher up in 
the disk photosphere. Large grains can exist much closer to the star than small grains since 
they cool more efficiently (grain settling in the inner rim with large grains dominating the 
mid-plane opacity may b e needed to explain near-IR SED and interferometry data on VVSer; 
Pontoppidan et al.l 120061 ) . A variation in grain size with disk scale height in the inner rim 
causes the rim to curve as shown in Figure [2l This curvature is on a much larger scale than 
the curvature produced in the IN05 models. Figures [7] and [9] show that the two models can 
be distinguished based on their visibility curves. 

In the dust segregation model that we studied in section 14.11 the effects of dust growth 
and sedimentation on rim structure were calculated assuming that the dust sublimation 
temperature is independent of density. But in reality, both dust evolution and the density 
dependence of dust evaporation temperatures (IN05 model) probably simultaneously play a 
role in determining rim structure. We have shown that in the absence of strong turbulence 
the rim shape is predominantly determined by grain growth and settling. Adding a density 
dependence to dust sublimation temperatures will lower rim height (see Fig. 2) and extend 
it further into the disk by just a few tenths of an AU: a small perturbation compared to the 
effect of dust segregation on rim shape. 

Calculations in this paper do not take into account effects of gas on dust rim geometry. 
Gas can affect the structure and location of the rim through accretion heating a nd by provid- 



i ng a so urce of opacity to stellar photons within the dust destruction radius. lAkeson et al. 



( l2005bl ) have shown that gas accretion heating might contribute significantly to the near-IR 
excess observed in some TTauri stars and play an impo r tant r ole in fixing the location of the 
dust destruction radius in these systems. lEisner et al.l (120061 ) argue that emission from gas 



within the inner rim is needed to explain spectrally dispersed near-IR interferometry data 
on Herbig Ae/Be systems. The effects of gas on dust rim structure and obse rved near-IR ex- 



cesse s in YSOs is not well understood. Some theoretical efforts (see for e.g. iMuzeroUe et al. 



20041 ) have been made in this direction but further progress would require detailed modeling 
of the relevant gas physics. 

The dust rim in YS Os is a unique test-bed for theories of grain evolution. A number 
of observational results Jp'Alessio et al.lboOll . bood : Ichiang et aPboOll : iRettig et al.]l2006h 
have found evidence for dust growth and settling in the outer regions (tens of AU scale) of 
YSO disks. In this work we show that long baseline near & mid-IR interferometry provide 
an opportunity for similar studies on the inner rim at sub-AU scale. Furthermore, if the 
gas within the inner rim were to emit significantly in the near-IR (gas emission is not well 
constrained and preliminary results from CHARA on AB Aur suggest that this might be 
important) , then the baselines of CHARA(~ 300m) are capable of partially resolving the 
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hot gas emission, at-le ast for a few Herbig systems. Combined with accretion line diagnostics 
( iMuzeroUe et al.ll2004l ). near-IR interferometry could help pin down the structure of the inner 
gas disk. 

In the future, we will compare SED and interferometry data with our model predictions 
for the inner rim. This will constrain dust growth, settling and turbulent mixing parameters 



that recent theore tical simulations (jJohansen fc Klahn l2005l : iFromang fc Papaloizoul 12006 



Turner et al.ll2006l ) have begun to predict. These simulations show that dust grains up to 
100 microns in size are well mixed with gas within one scale height. In the extreme case of gas 
and dust being well mixed at all scale heights, the rim shape will be controlled entirely by the 
density dependence of dust evaporation temperatures. Long-baseline near-IR interferometry 
can distinguish between various dust mixing scenarios and furthermore, it might even be 
possible to observe the time evolution of grain sedimentation processes by examining the 
shape of the inner rim as a function of stellar age. Detailed comparisons between model 
predictions and data will also provide insight into the role played by gas. 



6. Summary and Conclusions 



We have described a Monte Carlo approach to calculating the structure of the inner rim 
in protoplanetary disks. The obvious advantage of this method over an analytic approach 
is that it allows for a self-consistent determination of the entire disk structure. Any effects 
of spatially varying dust properties on disk structure (and the resulting effects on SEDs and 
visibilities) can also be easily studied. 

We have shown that dust growth and sedimentation may play a dominant role in deter- 
mining the structure of the inner rim. The difference in cooling properties between grains of 
different sizes, combined with sublimation and a vertical gradient in dust sizes in the disk, 
leads to inner rims that are extremely curved (over scales > 0.7 AU in cylindrical radius for 
an AO star, see Figure [2]). Synthetic images (Figures E] & ED and visibilities (Figures [7] & [HD 
for the rim models have been computed, highlighting differences between various theoretical 
models. 

The inner rim dominates near-IR emission from circumstellar disks. Studies in the 
near-IR thus isolate the rim from rest of the disk and are therefore extremely powerful 
tools for characterizing ri r n properties. In the last decade or so, near-IR in terferometers 
dMillan-Gabet et aLlEoogl . I2OO1I : lEisner et allboosl . booi lAkeson et aP boOSbl lah have been 
used to determine the sizes of inner rims in a number of Herbig Ae/Be and TTauri stars. 
This has helped in establishing a tight correlation between the radii of the rims and the 
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luminosities of the central stars (IMonnier fc Millan-Gabetll2002l : iMonnier et al.ll2005l ). The 



current generation of near-lR interferometers (CHARA, VLTI) are capable of taking the next 
natural step — determining the geometry of the inner rims. As shown in this paper and IN05, 
the geometry of the inner rim is tied to multiple dust properties such as grain evolution, 
grain sedimentation (which depends on the strength of turbulence), and dust sublimation 
temperatures (IN05). Once the inner rim shape is constrained, mid-IR interferometry (as 
seen in Figure [8], both the rim and the flared disk region emit in the mid IR) can be used to 
assess the geometry of the outer, fiared-disk component. 

A significant number of Herbig Ae/Be stars have multi-wavelength (near-IR, mid-IR 
and mm-interferometry, and spectroscopy over a large wavelength range) data available. In 
the near future we will try to generate disk models that simultaneously fit interferometry and 
spectral data over multiple wavelengths for a small sample of stars. This will be a stringent 
test for current disk models and should help alleviate model degeneracy that plagues single- 
wavelength studies. 



AT wishes to thank R. Kurosawa for parallelizing the image making modules in TORUS, 
A. Isella for providing data points for the IN05 curves in Figure [2] and N. Calvet for useful 
discussions. This project was partially supported by NASA grant 050283. This publication 
makes use of NASA's Astrophysics Data System Abstract Service. 



A. Appendix 

Here, we derive a general analytic formula for the rim shape of a dust segregation model 
[a case with specific numerical values for parameters, was discussed in section §4.2.21) . 

Let the densities of the two dust components in the circumstellar disk be given by 



PL 



arge 



^largePo_gas GXp 



1 



2 (Clarge^) 



(Al) 



Psmall — ^smallPo_gas exp 



1 



2 (Csmall^)' 



(A2) 



where piarge and psmaii are the densities of large and small grains respectively and Po_gas is 
the gas density in the mid-plane. Aiarge and Asmaw are normalization factors for the dust 
densities, z is the height from the disk mid-plane, h is the gas scale height and Ciarge and 
Csmaii are the scale heights for large and small grains as fractions of the gas scale height. 
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The ratio of Planck mean opacities at the dust evaporation temperature and the stellar 
effective temperature - e(z), then becomes 



v4iarge^Pl,_(Tevp) + Assail ^P,^,n (^cvp) exp 



1 -^'^ (Ismail ^largc) 
^largc ^small 



1 ^^(^small ^fargc) 



^large^small 



(A3) 



where Kp^^^^^ and Kp^^^^^ are Planck mean opacities for large and small grains at the spec- 
ified temperatures, Tevp is the dust evaporation temperature and T* is the stellar effective 
temperature. 



Neglecting local ac cretion heating, the dust destruction radius i?evp (see llsella et al. 



20061 : ICalvet et al.lll991[ ) is 



i?evp [AU] = 0.034 



m( 



1500 K 

\ cvp 



1 



e{z) 



1/2 



(A4) 



is luminosity of the star in solar units. B is a dimensionless diffuse heating term, charac- 
terizing the rim optical depth to its own thermal emission. 



The gas scale height h at Revp is estimated to be (see eqs. [7] and [8] from lChiang &: Goldreich 



19971) 



h [AU] = i?evp [AU] 



gas 



Mrr 



O 



cvp 



1/2 



(A5) 



^4 X lO^Ky \M^J V4.64 x 10-3 AU^ 
where Tgas is the temperature of the bulk of the gas (assumed to be isothermal) and is 
mass of the star in solar masses. 

Combining equations (A3), (A4), and (A5), the height of the evaporation front as a 
function of distance along the disk mid-plane from the star becomes 

Aiarge-f^Plargc(^*) " 'i^^large-f^Pi,,ge (^evp) 



Z[AU] =i?evp[AU] 



^^small-^Psmall(^cvp) — ^small-^Ps 



1/2 



X 







4 X io^k; 

where Y and Cdifr are dimensionless terms given by 



■evp 



4.64 X 10-3 AU 



1/2 



Y 



•l .71xlO-^°AU-^K-Xlp?;tp 



B 



(A6) 



(A7) 



Cdiff 



y-2 /■2 
SsmallSlarge 



Ismail 



- ^2 
^ large 



(A8) 
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Table 1. Basic properties of central star and the circumstellar disk 



Star Circumstellar Disk 



Mass 


2.5 Mo 


Surface Density 


S(r) = 2000(r/AU)-i-5 g cm-^ 


Teff 


10, 000 K 


Disk outer-Radius 


200 AU 


Luminosity 


47 Le 


Mass 


3.8 X 10-2 Mo 


Distance 


150 pc 
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(Optical depth to stellar photons) 



0.05 



3.98E-01 4.3gE-01 4.80E-01 5.21 E-01 5.62E-01 



0.2 
0.1 
-0.0 
-0.1 
-0.2 
-0.3 



Q. 1.0 1.2 1.4 1.6 1.8 2.0 



6.31 E-01 7.44E-01 8.57E-01 9.71 E-01 I.OBE-i-OO 




1.0 1.2 1.4 1.6 1.8 2.0 



(Optical depth to stellar photons) 



0.05 



6.31 E-01 6.88E-01 7.46E-01 8.03E-01 8.61 E-01 



J (b) 




1.2 1.4 1.6 1.8 



6.31E-01 8.14E-01 9.98E-01 I.I8E-1-OO 1.36E-I-00 



(d) 




1.0 1.2 1.4 1.6 1.8 2.0 



Distance (along mid-plane) from star[A.U.] 



Fig. 1. — a) Cross section of an inner rim from which dust has been stripped. Photons are 
propagated through this rim to determine the 'optically thin' grid-cell temperatures. The 
color scheme shows integrated r, measured along lines perpendicular to the disk mid-plane, 
b) Inner rim after the first dust growth step. Dust is grown in cells that are cooler than 
the sublimation temperature. The contours connect points with equal integrated tau. c) 
and d) depict stages further along in the dust growth scheme. The geometry of the rim is 
"sublimation converged" in (d) (see §2.21) . A hydrostatic equilibrium calculation (see §2.11) 
is then performed. The dust growth and hydrostatic equilibrium calculations are repeated 
until convergence is reached for the structure. The rim shapes depicted above are not the 
final converged solution. See Figure IB for the final rim shapes. 
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t[A.U] 


0.15 


o 


— IN05 values for density dependent dust sublimation model 
TORUS values for density dependent dust sublimation model 


/ = 




0.10 


I * 


— * TORUS values for dust-segregation model (this work) y 


^ / - 


■<1J 




— Analytic estimate for dust-segregation model ^ '' 


^/ ~ 


Rim 


0.05 




large ^ ^ - '\^^^(r^ dust 
grains ^ i^^js^S^^ — segregation 


^ small Jl 
^ grains 




0.00 










0.0 


0.5 Radius[A.U] 1 .0 


1.5 



Fig. 2. — The 'rim' is defined as the r = 1 surface (for A= 5500A), computed along radial 
lines from the central star. The figure shows the height of the inner rim above the disk mid- 
plane. The IN05 rim has been scaled at the ~ 8% level to match up with the TORUS rim. 
The dashed lines are an analytical estimate of the evaporation front for the dust segregation 
model (see gJ^D- 
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Fig. 3. — Near infrared emission (integrated between 1.25-7/im) from the inner rim as a 
function of inclination angle, plotted for small and large grains. The emission has been 
normalized to the stellar luminosity. 





Number of gas scale heights 



Number of gas scale heights 



Fig. 4. — a) The ratio of Planck mean opacities for dust at the sublimation temperature and 
at the stellar photospheric temperature for the dust segregation model, b) Effective dust 
size plotted against the number of gas scale heights (see §4.ip . 
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Fig. 5. — a & b show near and mid-IR SEDs of the star + rim system for the IN05 and dust 
segregation models, c & d show SEDs of the star + rim + the disk. The star is placed at 
150pc with stellar parameters described in Tabled! (a, c) system is face-on, (b, d) system is 
inclined 60*^ from face-on. 
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Fig. 6. — Synthetic 2.2/xm images for the different rim models discussed in the text. The 
panels on the left and right are IN05 rims computed for 1.2/im (large grain) and O.l/^m (small 
grain) silicate dust. The center panels are images for the dust segregation model. The star 
is placed at 150pc with the stellar parameters described in Table [TJ The star is unresolved 
at the image scale and is just one bright pixel at the center of the images. 




Fig. 7. — 2.2//m visibilities for the IN05 and dust segregation models. The panel on the left 
shows the visibilities for a facc-on disk and the right panel shows visibilities computed along 
the major and minor axis for an inclined disk. 
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Fig. 8. — Synthetic 10.7/xm images for the different rim models discussed in the text. The 
panels on the left and right are IN05 rims computed for 1.2/im (large grain) and O.l/xm (small 
grain) silicate dust. The center panels are images for the dust segregation model. The star 
is placed at 150pc with the stellar parameters described in Table [H The star is unresolved 
at the image scale and is just one bright pixel at the center. 
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Fig. 9. — 10.7//m visibilities for the IN05 and dust segregation models. The panel on the 
left shows the visibilities for a face-on disk and the right panel shows visibilities computed 
along the major and minor axis for an inclined disk. 



